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On the basis of requirements of unitarity and analyticity we analyze the real and imaginary parts 
of the pp’^ scattering amplitudes at recent ultrahigh energies, 1-100 TeV. The predictions for the 
region ^/s > 100 TeV and < 0.4 GeV^ are given supposing the black disk asymptotic regime. It 
turns out that the real part of the amplitude is concentrated in the impact parameter space at the 
border of the black disk. The interplay of hadron and Coulomb interactions is discussed in terms 
of the A-matrix function. The pp diffractive scattering cross section at 7 TeV is calculated with 
Coulomb interaction taken into account. 
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I. INTRODUCTION 

Ultrah^h energy proton-proton interaction data at 
LHC [H-Q and the cosmic ray data @ definitely tell 
us that at energies y/s ^5 — 50 TeV the asymptotic 
regime is switching on for diffractive scattering processes. 
These data together with those of ISR 0 demonstrate 
the steady growth of atot, o'ei and amei with the energy 
increase, a shrinkage of the diffractive cone in daei/dq^ 
and a relative supression of the real part of the scattering 
amplitude 

The growth of the cross sections atot, o'ei , Oinei and the 
diffractive scattering slope is consistent with the picture 
of fast moving hadrons as parton clouds with an increas¬ 
ing transverse size. The gluonic origin of the parton cloud 
explicates a slow growth of the diffractive scattering slope 
and a late start of the asymptotic regime: the effective 
mass of the soft gluon is notsmall being of the order of 
800- 1000 MeV l3,i]. 

The observed growth of total cross sections at preLHC 
energies 1,0 initiated studies of models with the su¬ 
percritical pomeron [HHIl- The discussion of the power 
growth of cross sections with energy actualized the prob¬ 
lem of s-channel unitarization of scattering amplitudes 
and the use the Glauber approach [l^. Taking into ac¬ 
count the s-channel rescatterings, the power-s growth of 
amplitudes is dampened to the Qn^ s)-type |16l - [la| . to the 
limits of the Froissart bound [Q . Still, let us emphasize 
that exceeding it does not violate the general constraints 
for the scattering amplitude [^ . 

The black disk picture appears to be a rather natu¬ 
ral mode for the ultrahigh energy corresponding to non¬ 
coherent parton interactions in hadron collisions. For 
the black disk scenario the profile function at i/s > 100 
TeV gets frozen inside the disk area, T{b) ~ 1 at 
b < Rbiack disk , while the increasing radius of the black 
disk, Rbiack disk , determines the total, elastic and inelas¬ 
tic cross sections: atot ^ ‘^T^Ruack disk^ ^ei ^ T^Rhack disk 
and atnei — '^Rliack disk ' black disk mode was in¬ 


tensely discussed in the last decade, see, for example, 
[2ll - l^ and references therein. 

In the present paper we underline that at asymptotic 
energies the imaginary part of the amplitude Ag (q2,lns) 
turns into a generating function for the real part, 
Ag^(q^,lns), due to unitarity and analyticity require¬ 
ments. In Section 2 the real parts of the hadronic scat¬ 
tering amplitude are calculated for a set of energies, 
•y/s = 1,10,10^,..., 10® TeV, and profile functions, Ta(&) 
and T^(b), are presented. In Section 3 we discuss a com¬ 
bined action of the Coulomb and hadronic interactions 
for the diffractive scattering region. If here the eikonal 
approach works, the straightforward way to take into ac¬ 
count the interplay of these interactions is the use of the 
AT-matrix function technique, thus keeping valid the uni¬ 
tarity condition. We present the corresponding formulae 
and calculate the pp diffractive scattering cross section at 
7 TeV with the Coulomb interaction taken into account. 


II. REAL PART OF THE HADRONIC 
SCATTERING AMPLITUDE 

The high energy scattering ampitude is dominantly 
imaginary, the real part of the amplitude is a next-to- 
leading term. So, in the pre-asymptotic region we should 
include into consideration of the scattering amplitude 
both the imaginary and real parts. The amplitude reads: 

^(q^e) = J d%c^'^T{b,o, 

T{b,0 = TQ{b,0-^Mb,0, ( 1 ) 

where ^ = Ins, & = |b|. For the profile function we write: 

T{b,^) = l- r]{b, 0 exp {2i5{b, ^)) = , (2) 

1 - iK{b,Q 

presenting it in terms of the phase shift 6(b,^) and in¬ 
elasticity parameter r]{b,^), or the AT-matrix function 

Kib,0- 
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FIG. 1: Profile functions for imaginary and real parts of the pp scattering amplitude, T^{b, and T^{b, ^), at a set of energies 
^/s = 1,10,10^,..., 10® TeV; for ^/s > 100 TeV the black disk mode is suggested. 




FIG. 2: Imaginary and real parts of the pp scattering amplitude, and ^), at energies y/s=l (solid),10 (dashed), 

10^ (dot-dashed), 10® (dotted) TeV; for >>10 TeV the black disk mode is suggested. 


If the imaginary part of the amplitude dominates, the 
next-to-leading terms include real part due to analytic- 
ity requirement. Namely, we should take into account 
contributions both of the s-channel and the w-channel. 
Suggesting atot{pp) = crtotipp) at s -)> oo we write: 


1 

2 1' 


A(q^lns)-b A(q^ln(-s))J ~ 
9A(q^,lns) —in 


~ A(q , Ins) -b 


~ AQ(q^,lns) -b 


i9(lns) 2 
9AQ(q^,lns) —in 
9(lns) 2 


(3) 


where ln(—s) = Ins — in. It means the amplitude 
^9(q^ Ins) is the generating function for Ajf(q^,Ins). 
Analogously we write for the profile function: 


T3j(5,lns) 


n dTQ{b, Ins) 
2 (9(lns) 


(4) 


The usual notation reads Tjf(&, ^)/Ta(&, ^) = p{b,^), 


therefore the total and elastic cross sections are written 
as: 


Tib,0 = (l + ^p)^Q(6,o, 

atot = 2 J d^bT^ib,^), 

47r^ = (1 + P^)^|(q^)- (5) 

Taking into account that is small, p^ ~ 0.01, one can 
approximate: 


bl^(q^0 



( 6 ) 


that make possible direct calculations of the real part 
of the scattering amplitude, A 3 j(q^,^), on the basis of 
the energy dependence of the diffractive scattering cross 
section. 
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FIG. 3: Imaginary and real parts of the Jf-matrix function, and Ksii{b,^), at energies y/s = 1,10,10^,..., 10® TeV. 


A. Calculation of the real part of the scattering 
amplitude 

Using Eq. dO]) we calculate and, correspond¬ 

ingly, the profile functions T^{b, and TQ{b, £), see Figs, 
mu The figure [1^ for TQ{b, £,) is taken from where 
a comparison with data is given as well. Calculations at 
> 100 TeV are done suggesting the black disk mode 
as a realization of the asymptotic regime. 

An advantage of the AT-matrix function technique is 
the separation of the two-particle rescattering states 
which turn out to be mass-on-shell for leading terms of 
the amplitude [l^. Following Eq. , we write: 

= ■ (7) 

The functions K^{b,^) and KQ{b,^) for i/s = 
1,10,10^,..., 10® TeV are shown in Fig. 


B. Eikonal approach for scattering amplitude at 
ultrahigh energies and the Feynman diagram 
technique 

For the pp-scattering amplitude App_>pp(ppi„ —>■ ppout) 
the reproducing integral reads: 


{Po,P±,Pz) = (P + w^/2p,0,p) and p 2 = {p + 
m^/2p, 0,—p). For intermediate and final state protons 
we have: 


k;i'_L + k2'_L = 0, ki^ -1- k2_L = 0. 

92' = iP2 - k2'f -k|_L, 

92'2 = (*2' - k^f -(k2i - \^2'±f , (9) 


where ki,k 2 refer to momenta of outgoing protons. 

At ultrahigh energies the AT-matrix function is dom¬ 
inantly imaginary for the black disk and resonant disk 
modes . That means the dominance of the mass- 

on-shell contribution in the loop diagrams. For the 
rescattering diagrams this is realized in the replacement: 


(m^ — ky — i0)(m‘^ — A:|, — fO) 
—27r^S(m^ — ky)d{m^ — A:|,) = 

= —27r^d (jt[t^k[7^ — {w? -\- kpj_)^ 

xd — {w? -\- k2/j^)^ , 




( 10 ) 


App—>pp (pPm ^ PPout') — ATpp—^pp (ppin ^ PPout) 

7fpp->pp (1^2' —» PPout) 

(m^ — ky — i0){m? — — 10) ’ 

to be definite we consider proton-proton scattering. Here 
Kpp^pp is the block without two-particle states thus be¬ 
ing up to factor the AT-matrix function, indices (1', 2') 
refer to protons in the intermediate state. 

Let us consider the scattering amplitude in the cm- 
system where we write for the initial protons: pi = 


where fcA) = ko + kz, k^ ^ = ko — kz- Then the right- 
hand side of Eq. m reads: 

A2-!.2(k2_L)C) = -^2-).2 (k2j_,‘f) -b (H) 

+ J (27r^)^ A2^2(k|j_,^)—A:2^2 ((k2'_L - k2_L)^,C) , 

here Ar2^2(k5_, C)/(4s) = Ar(k^,^) is the AT-matrix func¬ 
tion in momentum representation. 
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[3=;ib [3=;ib 

FIG. 4: Pure coulombic interaction in pp collision at different A: the function K'^ (b) in hadronic region of the impact 
parameters, 6 < 25 GeV”^, and coulombic region, (3 = Xb > 1. 


C. Impact parameter presentation 

The Fourier transform gives the iF-matrix function in 
the impact parameter space: 

■^^ 2^2 = J (fbexp{{kh)K{b,^), 

^^ 2^2 (ki,C) = y fi^6exp(fkb)a(6,^), (12) 


A. Interplay of hadronic and Coulomb interactions 
in the A-matrix function techniqne 

We consider two types of scattering amplitudes and 
corresponding profile functions: the amplitude with com¬ 
bined interaction taken into account, ,^) and 

^), and that with switched-off the Coulomb 
interaction, following to © we use for them notation 
Al(q^,^) and T{b,^). For the combined interaction pro¬ 
file function we write: 


Equation (HU in the impact parameter space is written 
as: 


a{b,0 = lK{b,0 + a{b,0^K{b,0■ ( 13 ) 


r^+^(&,e) 


-2fA:C+^(&,^) 
l-iK(^+^(b,0 
-2* {K^'{b)+K{b,0) 
l-i{KC{b)+K{b,0) ’ 


Thus, we have the formula of the eikonal approach: and the Coulomb interaction is written as: 


(15) 


a(6,C) 


iK{b,0 

1 - iK{b,C)' 


(14) 


The function K{b,^) depends on the energy and realizes 
effectively the interaction which manifests itself in the 
shrinking of diffractive cones with the energy increase. 


III. DIFFRACTIVE SCATTERING 
AMPLITUDE AT ULTRAHIGH ENERGIES AND 
COULOMB INTERACTION 


\Cf2\ /„2\ 47ra 2\ 


A^(q^)=±*/i(q^)^, 

-2iK^{b )=±* y 


q- + A2 
cPq 


/2(q^) 


(16) 


(27r)^ 




Here a = 1/137; the upper/lower signs refer to the 
same/opposite charges of the colliding particles. The cut¬ 
ting parameter A, which removes infrared divergency, can 
be put to zero in the final result for A^+^(q^,^). Col¬ 
liding hadron form factors, /i(q^) and / 2 (q^), guarantee 
the convergency of the integrals at q^ —oo; for the pp^- 
collisions we use: 


The interplay of hadronic and Coulomb interactions 
was studied in a set of papers, see [sol - ls^ and refer¬ 
ences therin. At ultrahigh energies and small q^, where 
the eikonal works, the straightforward way to take into 
account the combined action of hadronic and Coulomb 
interactions {H + C) is to use the technique of the K- 
matrix function. Here we illustrate this way by corre¬ 
sponding calculations of {b,^) and profile function 

r^+'^(&,e). 


/i(q^) = / 2 (q^) 


1 

(1+071^)^ 


(17) 


B. Numerical calculations 

In Fig. 0] we show K^{b) for A = 0.1 CeV and 
0.01 GeV. The inclusion of the Coulomb interaction into 
consideration of hadron diffractive scattering does not 
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FIG. 5: Combined hadronic and coulombic interaction for the pp scattering: Real parts of the iF-matrix functions at different 
A at a set of energies ^/s = 1,10,10^,10® TeV; for ^/s > 100 TeV the black disk mode is suggested. 




FIG. 6: Gombined hadronic and coulombic interaction for the pp scattering: Imaginary parts of the proHle functions at 
different A at a set of energies y/s = 1,10,10^,..., 10® TeV; for -^s > 100 TeV the black disk mode is suggested. 


change the imaginary part of the iC-matrix function, 
{b,^) = . The real parts of the iC-matrix 

functions {b,^) = K^{b,^) + K^{b) for different A 

are shown in Fig. O 

Imaginary and real parts of the profile functions, 


f) for the hadronic region, 5 < 25 GeV“^, are 
shown in Fig. [51 considerable perturbations are seen in 
the real part. 

With these A’s we calculate at ^/s = 7 TeV the 
profile function {b,^LHc) and the corresponding 
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FIG. 7: Diffractive scattering cross section for pp at 7 TeV [I[ versus description with interplay of the Coulomb interaction 
(we use here A=0.01 GeV) and hadronic one (the real part of the hadronic amplitude is taken into account): figures (a,b) refer 
to version (1) for determination of the hadronic amplitude, figures (c,d) to version (2); solid curves refer to pp, dashed ones to 

pp. 


amplitude The determination of the 

hadronic amplitude, ^^.(q^, ^lhc), is performed in terms 
of two versions: 

1) with a direct application of the approximation ([5]) to 
the TOTEM data Q, 

2) using the the results of the Dakhno-Nikonov model 

HilSl. 

The description of the data for in terms 

of these two versions is shown in Fig. [7j here Figs. [7^,b 
refer to the version (1) and Figs. [71:,d correspond to the 
version (2). 

The Dakhno-Nikonov model gives a somewhat worse 
description of the gg_ 

ing Eq. ®. This is not surprising because the model 
descibes the data in a broad energy interval, 0.5-50 TeV 
, and the model parameters are responsible for a com¬ 
plete set of the data. 

The specificity of the if-matrix function treating the 
amplitudes with the Coulomb interaction is the use of 
the determination: (h) = tg6^{s,b). Another deter¬ 

mination was applied in [^: K^{b) = S^{s,b). 


IV. CONCLUSION 

On the basis of requirements of analyticity we calculate 
the leading terms of the real part of the pp-scattering am¬ 
plitude for diffractive interactions at ultrahigh energies, 
^/s > 1 TeV, and small momenta transferred, < 0.4 
GeV^. We do not include into consideration the region 
with larger q^: in the region of larger momenta trans¬ 
ferred mechanisms with conventional Pomeron as well as 
short-range non-Pomeron interactions are possible (for 
example, see Refs, but here we concentrate our 

attention on peripheral interactions. 

In the region of the diffractive scattering cone the imag¬ 
inary part of the amplitude prevails over the real part 
Agif(q^,^)/Acj(q^,^ 1/^. The unitarity and analyt¬ 
icity requirements give unambiguously the leading term 
of the real part. We calculate Ag^(q^,^) at q^ < 0.4 
GeV^ supposing for energies ^/s > 100 TeV the black 
disk mode. 

Presently we have a number of papers devoted to the 
asymptotic behavior of diffractive amplitudes at ultra- 
high energies, see for example |4Ill - l4^ and references 
therein; in these papers, however, the analyticity con¬ 
dition related to the u-channel is disregarded. 

The interplay of the hadronic and Goulomb interac- 
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tions at very small q_L^ is discussed in terms of the K- 
matrix function. The specificity of the scattering ampli¬ 
tude at ultrahigh energy is the dominance of the mass- 
on-shall contributions in intermediate rescattering states 
that results in the mass-on-shell origin of the ilT-matrix 
functions. Such ilT-matrix functions allow to incorporate 
the Coulomb interaction terms into the scattering am¬ 
plitude straightforwardly, keeping te unitarity condition. 
We present corresponding formulae and perform calcula¬ 


tions for the black disk mode. 
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